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Since 1908, when Mie reported analytical expressions for the fields scattered by a spherical particle 
upon incidence of an electromagnetic plane-wave, generalizing his analysis to the case of an arbitrary 
incident wave has proved elusive. This is due to the presence of certain radially-dependent terms in the 
equation for the beam-shape coefficients of the expansion of the electromagnetic fields in terms of vector 
spherical wave functions. Here we show for the first time how these terms can be canceled out, allowing 
analytical expressions for the beam shape coefficients to be found for a completely arbitrary incident field. 
We give several examples of how this new method, which is well suited to numerical calculation, can be 
used. Analytical expressions are found for Bessel beams and the modes of rectangular and cylindrical 
metallic waveguides. The results are highly relevant for speeding up calculation of the radiation forces 
acting on spherical particles placed in an arbitrary electromagnetic field, such as in optical tweezers. 

PACS numbers: 03.50.De,41.20.-q,42.25.-p,42.25.Bs 



Gustav Mie, in his celebrated 1908 paper fH], used the 
vector spherical wave function (VSWF), or partial wave 
expansion (PWE), of a linear polarized plane-wave to gen- 
eralize scattering theories to spherical particles of any size, 
from geometrical optics to the Rayleigh regime, and thus 
was able to clarify many phenomena, for example in atmo- 
spheric physics. He obtained analytical expressions for the 
expansion coefficients based on special mathematical iden- 
tities related to a plane-wave. This beam expansion was 
necessary for applying boundary conditions at a spherical 
interface. Since then, with the arrival of lasers and optical 
waveguides, the diversity and complexity of possible inci- 
dent fields has become enormous so that the restriction to 
an incident plane-wave has become unrealistic. 

Different experiments, ranging from particle levitation 
and trapping [z,^, to the ultrahigh-Q microcavities used 
in cavity QED experiments [4, 5], use different beams. For 
example, very high numerical aperture beams are used in 
optical tweezers and confocal microscopy IB-Si) evanes- 
cent fields in near-field microscopy [9, 10], and the waveg- 
uide modes of a fiber taper are employed to couple light 
to the whispering gallery modes of spherical microcavi- 
ties iflUl. Optical forces, absorption, Raman scattering and 
fluorescence can be greatly enhanced inside spherical mi- 
crocavities at Mie resonances lfl2l - [l5ll . Laguerre-Gaussian, 
Hermite-Gaussian and Bessel beams [16, 17], and the in- 
ternal electromagnetic field of hollow core photonic crys- 
tal fibers fTS, T9], are used to trap and transport particles. 
The understanding of all these phenomena requires a pre- 
cise knowledge of the VSWF coefficients of the incident 
beams. A generalized Lorenz-Mie theory was developed to 
handle the many variants of beams beyond classical plane- 
waves, and the expansion coefficients in these cases are 
known as beam shape coefficients (BSC) lEolElll . More- 



over, because the VSWFs form an orthogonal complete ba- 
sis, they can be used to study scattering and forces ll22ll 
on non-spherical particles, and are the starting point of the 
powerful T- matrix methods f23\. 

The calculation of BSCs for an arbitrary beam has al- 
ways been a complicated task, requiring significant effort. 
Furthermore, there is a fundamental problem with these 
calculations: an expansion of any function in some basis 
is complete only when the expansion coefficients can be 
written in terms of scalar products, or integrals, with de- 
fined numerical values. This task has actually never been 
accomplished for the VSWFs of an arbitrary beam because 
the integral over the solid angle does not explicitly elimi- 
nate the radial dependence, at least up until now. So far as 
we are aware, the current literature lacks any mathematical 
proof that this radial function, which appears after integra- 
tion over all solid angles for any type of beam that satisfies 
Maxwell's equations, can exactly cancel out the spherical 
Bessel function that appears on the other side of the BSC 
equation. If this is not true, then the BSC could not be a 
constant independent of the radial coordinate - as required 
for a successful expansion. 

This non-radial dependence of the BSC has been proven 
only for the case of plane-waves and for a high numeri- 
cal aperture focused Gaussian beam fl^. Working with an 
electromagnetic mode inside a hollow cylindrical waveg- 
uide, we also have been able to obtain analytical expres- 
sions for constant BSCs that depend only on the position 
of the reference frame. This raises the fundamental ques- 
tion, whether it would be possible to prove that the spheri- 
cal Bessel function would naturally emerge from the solid 
angle integral for any type of electromagnetic field. The 
purpose of this letter is to show, we believe for the first 
time, that this is indeed possible. The implications of this 



2 



result for computational light scattering is very noteworthy. 
We show how the new method can be used to calculate the 
BSCs for plane-waves, cylindrical and rectangular waveg- 
uide modes and Bessel beams. 

The dimensionless BSCs GL^/^*^ for an incident field 



So we define the k-only dependent fields 
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E = E(r), H = H(r) are defined in the equations |,25|] Imposing this restriction one finally obtains that 
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where En is an electric field dimension constant, /cN 
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f), jp{kr) are spherical 
LYpg{r)/,/p{p+l) are the 
/Jjje, k = Uy/JIe and L = 



iV X Mpg, Mpg = jp{kr):Kpg 
Bessel functions, X(f) 
spherical harmonics, Z 
—ir X V. Throughout this paper we use the convention 
that the terms inside [] are parts of separate equations, i.e., 
(d} contains two equations, the first relating E to Gp^ and 
Gjf and the second, ZH to G™ and -Gjf . 

The usual procedure for obtaining the BSCs involves 
multiply both sides of ([T]l by X*,^,, take scalar products 
with the fields and integrate over the solid angle $7. Due 
to the orthogonality properties of the vector spherical har- 
monics 125112^ . one can easily show that 
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where 17 is the solid angle with respect to an arbitray 
origin not related to any particular point of the incident 
beam. Equation ^ does not yield explicit expressions 
for the BSCs because the LHS still contains the radially- 
dependent spherical Bessel function. Our goal is to extract 
this function from the RHS and cancel it out with the one 
on the LHS, for any general incident electromagnetic field. 
To accomplish this we use the Fourier transform of the 
fields 
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By this definition, one can show that J^{Li/)(r)} = 
C^{]s.) and that the angular momentum operator in re- 
ciprocal k-space C has the same form as in real r-space 
(is Hermitian) and is given hy C = —ik x where 
V/c = {dk^,dk , dk^)- Using this property and the expan- 

sione"^'-- = 47rEz=o^'jK^'OEL-/^™(f)>^™(k) 
and making A4pq(k) = jp{kr)Xpg{k), Xpq(k) = 
CYpq{k)/ ^Jp{p + 1), we obtain 
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Now, only Fourier transforms of the form fc^^(k') = 

b{k' — A;)^k(k ) will represent a field i^(r) that satisfies 
the wave equation V^i^ + k-^F = in three dimensions. 
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To calculate the coefficients placed at arbitrary position 
Tq independent of the coordinate system of the fields one 
can use the translation property of Fourier transform [26]. 
Since this form is free of any radially-dependent function, 
our goal has been achieved. 

To calculate these Fourier transforms numerically it 
could be convenient to use Laplace series ll2ql . also 
known as spherical harmonic transforms ||27[|. Us- 
ing the Laplace series expansion [£{k), Z'H{k)] = 
Ep',„' [ep',q', V,9']^p*,g'(k ) one obtains 
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The e;/ „j//h;/ „j/ coefficients can be calculated by sev- 
eral algorithms freely available in internet and the matrix 
(p', q'\C\p^ q) is shown in most quantum mechanics books. 

Plane wave: The fields of a linearly polarized plane 
wave are given by [E(r),ZH(r)] = E^le^k x e]e'''^''", 
where e is the unit polarization vector, therefore the angu- 
lar Fourier transform is easily calculated, one obtains from 
equation Q that 
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For the special case k = z and for a circularly polarized 
wave with e = x it iy we have 
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where 5„ ,„ is the Kronecker delta function. This is the 
result shown in Jackson's book I25I1 . 

General hollow waveguide mode: The TM and TE 

modes of a hollow, cylindrical waveguide of arbitrary 
cross-sectional shape are given as function of g{r) = 
gp{p)e'^'^'^ , where gp{p) is the scalar solution of the trans- 
verse wave equation + dy + 7^)5p(p) = satisfying 
the boundary conditions at the waveguide surfaces 12511 . 
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where is the wavevector in the z direction, k'^ = k'^+j'^, 
7 is the transverse wavevector. 

It should be useful to introduce the cylindrical coordinate 
system in the r-space and k-space. In the r-space we have 
p = XX + yy, = p ■ p and p = p/p. We have also 
cf) = — yx+xy and0 = cf)/p. In the same way, ink-space 
we change p — )• 7 and (f) ^ (j)^, obtaining an equivalent 
system of coordinates. So, the Fourier transforms of fields 
([Tol l and (fTTb are given by 
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Now, similar to the previous case, only Gp(k') = 
G^{(p'i^)5{'j' — 7)/7 can represent in the Fourier domain 
a scalar function gp{p) that satisfies the transverse wave 
equation (d^ +dy+ 7^)5p(p) = 0. So, we have G'(k') = 
G,{cP',)6{y-^)6{k'^-k,)/i = G,{(t>',)5{k' -k)6{e',- 
9k)/{sm0'f^k''^). Applying these results to and ([T3l l 
and using ^ we obtain 
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is an Associated Legendre Polynomial 

and cos 9k = k^/k, P^'{x) is the derivative of Pj'(x), 
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Rectangular hollow metallic waveguide mode: In 

this case there are two scalar functions gp, given by 
= cos(A:rx') cos(A:,"y') and gj''{p) = 
sin(/c™x') sin(/c^y'), where /c™ = rmr/a and kj^ = 
nir/b, (/)™" = arctan(/c^Y/c™) = arctan(na/m6) and 
7mn = T^\Jvn? jcfi + IT? /V^, fTi, u integers and the ori- 
gin at the lower left corner of the waveguide [251]. If 
the origin is placed at (xq, yo, 0) then x' = x + Xq and 
y' = y + yo- The Fourier transforms of these fields 
can be easily calculated, and the result is a sum of four 
Dirac delta functions a the points {k^,k^), {—k"\k^), 
(/c™, -A;^), {-kl\ -k]^). These delta functions can be 
written in cylindrical coordinates to give us the function 
G^{4>'^.). The integral over (j)'^. given by in ( [TTl l - where 
the plus sign is for the TE-waveguide mode and the minus 
sign is for the TM-waveguide mode - is 
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The BSCs can now be found substituting (fTSl l in (fT4] | and 

csi. 

Cylindrical hollow metallic waveguide mode: The 

scalar solution for the electromagnetic fields in terms 
of cylindrical coordinates with the origin on axis 
for the metallic waveguide are given by gp{p) = 
Jmi'ymnP)^^^™'^ , whcrc Jm{x) are order m Bessel func- 
tions, J^(x) is the derivative of J„(x), 7,„_„ = Xm,n/R 
or 7m, n = x'm nl P^ ^i^h R being the radius of the cylin- 
der, Xm.n the n-th root of Jm{x) for the TM mode and 
X'„ „ being the n-th root of J^(x) for the TE mode [25]. 
From now on we define the functions 
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and remember the addition theorem (l^ to express the 
scalar function in terms of a new coordinate system r = 
r' + Fq, written as a convolution 

00 

'0„(k;r)= ^ V'„_j(k;ro)V'j(k;r'). (21) 

j = -oo 

The Fourier transform of ^/^„i(k;r) can be written as 
^,„(k;k )6{k' - k)/k'^ and this way 
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The BSC'S can now be found substituting ( |23] | in ([T4l l and 
( [15]) . The on-axis case can be obtained by setting Pp = 0, 
which implies that Jm^q^jmnPo) = ^m.±q, and therefore 

G« = (2^)3/2(-i)±™5„,±,. (24) 

As we can see the sum over g in ([Tjl will disappear. 

Bessel Beams: We have calculated BSCs for two kind 
Bessel beams electromagnetic fields that obey Maxwell 
equations ll29ll . They are derived using vector potential and 
Lorentz gauge ^30] and are given as function of already 
defined in ^ and (jZB- 
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The Fourier transform of these fields are given by 
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The k' component is obviously null by orthogonality with 
angular momentum. Using again the addition theorem (l2T] | 
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and making c^f = \/p{p + 1) — q{q =b 1) one can show 
that the coefficients are given by ( [29l ) and dSOl l. 

In conclusion, we have shown that radially-independent 
amplitudes (the BSCs) of a complete set of vector spheri- 
cal wave functions can be calculated explicitly for an arbi- 
trary electromagnetic field. We have shown how this result 
can be used to determine the BSCs for several beam-types 
commonly employed in photonics, although of course the 
method is not restricted to applications within the field of 

I 



optics. This new-found ability to evaluate the BSCs of 
the VSWFs analytically makes it much easier to explore 
rapidly the influence of experimental parameters in prac- 
tical field scattering and optical tweezer systems. With 
this analytical breakthrough, the long-standing problem of 
evaluating the BSCs for an arbitrary field has been solved 
and the non-radial dependence of the BSCs proven, allow- 
ing one to avoid unnecessary approximations in the numer- 
ical evaluation of these quantities. 
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